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摘要 :通过 改进 的 光 能 利用 率 CASA 模型 估算 2001 一 2010 年 间 下 勒 河流 域 陆 地 生态 系统 的 净 第 一 性 生产 力 ( NPP ) ,采用 线性 
趋势 分 析 、 变 异 系数 和 Hurst 指数 等 方法 ,分 析 了 NPP 的 时 空 变化 特征 及 其 与 气候 因子 的 相关 性 。 结 果 表 明 : 中 蚊 勒 河流 域 
NPP 的 空间 分 布 具 有 明显 差异 ,空间 上 呈现 西北 低 .东南 高 的 趋势 ,日 具 有 较 明 显 的 经 向 “条 带 ” 分 布 特征 ,2001 一 2010 年 , NPP 
平均 值 为 102.26 gC ma'。@2001 一 2010 年 , 玻 勒 河流 域 NPP 总 体 呈 增长 趋势 ,年 际 波动 较 大 ,NPP 明显 增加 的 区 域 占 总 面 
积 25.15% ,明显 减 小 的 区 域 约 占 11.93%。@Hurst 指数 分 析 显 示 , 玻 勒 河流 域 NPP 变化 的 同 向 特征 强 于 反 向 特征 ,其 中 持续 改 
善 地 区 占 总 面积 的 78.3% ,21.7% 的 区 域 将 由 改善 转 为 退化 。 轴 在 年 尺度 上 ,降水 是 植被 NPP 变化 的 主要 影响 因素 ,NPP 与 降 
水 呈 弱 的 正 相关 关系 ,与 温度 相关 性 不 显著 ;在 月 尺度 上 ,温度 是 NPP 变化 的 主要 控制 因子 。 玻 勒 河流 域 NPP 对 气候 因子 不 存 
在 明显 的 时 淆 和 累积 效应 。 
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Abstract: Vegetation net primary productivity ( NPP ) is the solar energy fixed by vegetation during photosynthesis, which 
can be described as the net increase in biomass per unit time per unit area. The NPP, which plays an important part in the 
terrestrial carbon cycle, directly reflects the production capacity of terrestrial vegetation under natural conditions, and 
determines the ability of ecosystems to store carbon. The Shule River basin, located in northwest Gansu Province, is an 
important water source and ecological division within the Hexi Corridor. Increasing human populations and regional economic 
development have resulted in severe environmental problems within the study area, including natural disasters , vegetation 
degradation, and soil erosion. The study of NPP in the Shule River basin can provide a scientific basis not only for the 
prevention and control of natural disasters but also for the improvement of the ecological environment. This is important for 
local sustainable ecological development in the context of global climate change. The NPP in the Shule River basin was 
calculated using Moderate Resolution Imaging Spectroradiometer ( MODIS) images, meteorological data, and vegetation 


type maps, along with an improved Carnegie-Ames-Stanford Approach ( CASA) NPP model. The spatio-temporal 
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characteristics of NPP in the Shule River basin were analyzed for the period 2001 一 2010 using linear trend analysis the 
coefficient of variation, and the Hurst index. The relationship between NPP and climatic factors was investigated. The 
results showed that: (1) NPP in the Shule River basin had considerable regional differences. NPP gradually decreased from 
southeast to northwest in the Shule River basin, and presented a noticeable linear pattern in a longitudinal direction. The 


2 


inter-annual variability showed overall growth from 2001 to 2010. The annual NPP value ranged from 99.97 to 127.2 eCm 
a , with the minimum value occurring in 2001 and the maximum in 2010. The average annual NPP was 125.08 Cm a’ 
over the period 2001 一 2010; (2) however, within this overall increase, the annual NPP across the Shule River basin 
varied from 2001 to 2010. Of the total area of the basin, 25.15% showed a significant increase in NPP, whereas 11.93% 
showed a significant decrease; (3) the Hurst index analysis indicated that the majority of changes in NPP followed the same 
direction; 78.3% of the entire study area was predicted to follow a continuously increasing trend, while 21.7% of the total 
area will decrease in future; and (4) there was a weak relationship between NPP and annual precipitation, and no 
significant relationship between NPP and average temperature at an annual scale, indicating that precipitation was the key 
influencing factor. However, at a monthly scale, temperature was the dominant climatic factor affecting NPP in the study 
region. Analysis of the relationships between NPP, precipitation, and temperature showed that correlations at an annual 
scale were lower than those at a monthly scale. Precipitation and temperature exerted a variety of influences on different 
parts of vegetation NPP. The relationship between NPP and climatic factors suggested that NPP displayed no obvious time 


lag or cumulative effect due to climatic factors. 
Key Words: net primary productivity ( NPP); CASA model; temporal and spatial variation; MODIS ; Shule River basin 


植被 净 初 级 生产 力 (NPP) 是 单位 面积 和 单位 时 间 内 绿色 植物 从 光合 作用 所 获取 的 有 机 物质 总 量 中 去 除 
自 养 呼吸 以 后 的 干 物 质量 '"| 。NPP 能 够 直观 地 反映 出 地 表 植 被 在 自然 环境 状况 中 的 生产 能 力 , 常 常用 作 判 
定 碳 源 / 碳 汇 以 及 区 域 生 态 支持 能 力 的 重要 因子 '。 此 外 ,NPP 的 变化 对 温度 .降水 等 气候 因素 较为 敏感 ,而 
这 些 因素 与 植物 代谢 密切 相关 ,因此 NPP 的 变化 能 够 很 好 地 表征 全 球 气候 变化 。 掌 握 区 域 NPP 年 际 间 的 
定量 时 空 演 变 特征 趋势 ,对 于 理解 全 球 气候 变化 对 陆地 生态 系统 植被 变化 过 程 的 作用 机 制 ,科学 评价 区 域 生 
态 系统 的 调节 过 程 以 及 估算 陆地 碳 源 / 汇 具有 重要 的 理论 和 实践 意义 。 自 Ebermayer'”] 在 1876 年 最 时 开始 研 
究 NPP 以 来 ,NPP 研究 经 历 了 从 传统 的 站 点 实测 法 到 大 量 的 统计 回归 法 再 到 多 源 数 据 的 机 理 模型 法 3 个 阶 
段 。 在 区 域 尺 度 上 ,全 空间 覆盖 测量 NPP 显然 具有 极 大 难度 。 随 着 遥感 技术 的 发 展 , 利 用 模型 结合 遥感 数据 
进行 植被 NPP 的 估算 及 其 分 布 研究 逐渐 成 为 主流 的 研究 方法 。 目 前 国内 外 提出 用 于 估算 NPP 的 模型 多 达 
20 余 种 :5 ,大 体 上 可 归纳 为 3 类 :9 :气候 相关 模型 (统计 模型 ) . 光 能 利用 率 模 型 (参数 模型 ) 和 生态 系统 过 程 
模型 ( 机理 模型 ) 。Carnegie-Ames-Stanford Approach (CASA ) 模型 是 属于 基于 光 能 利用 率 机 理 的 一 种 过 程 模 
型 , 它 采 用 光 能 利用 率 与 光合 有 效 辐 射 两 个 因子 并 引入 气温 、 降 水 等 光合 作用 的 胁迫 参数 来 计算 NPP ,在 区 
域 和 全 球 陆 地 植被 NPP 监测 中 被 广泛 应 用 1。 在 区 域 尺度 上 ,由 于 植被 NPP 变化 受 植被 自身 生理 特征 差异 
的 影响 ,对 气候 、 地 形 、 土 壤 以 及 人 类 活动 等 因子 变化 的 敏感 性 和 响应 特征 存在 较 大 差别 , NPP 时空 演 变 的 驱 
动机 制 是 研究 的 难点 之 一 。 在 国内 ,学 者 们 针对 NPP 的 测算 方法 .时空 分 异 格局 驱动 因子 等 不 同 视角 开展 
了 研究 “” 。 研 究 区 域 多 集中 在 全 国 以 及 东北 .华北 等 大 区 域 或 省 域 尺度 ,而 针对 作为 一 个 较 完整 生态 系 
统 的 流域 尺度 的 研究 则 相对 较 少 。 

玻 勒 河 流域 深 居 内 陆 , 气 候 干 旱 , 生 态 环境 敏感 脆弱。 近年 来 ,受气 候 变化 和 人 类 活动 加 剧 的 影响 ,流域 
出 现 了 草场 退化 .湿地 蓉 缩 .土地 沙漠 化 次 生 盐 渍 化 .地 下 水 位 下 降 等 一 系列 生态 环境 问题 '… ,并 与 河西 走 
廊 其 他 内 陆 河 下 游 一 样 ,成 为 我 国 沙 尘 源 区 之 一 。 区 内 著名 的 风景 名 胜 月 牙 泉 由 于 水 位 下 降 严重 , 目前 靠 
“打点 滴 ? 来 维系 。 如 何 实现 流域 自然 -经 济 -社会 的 协调 和 可 持续 发 展 , 是 值得 探讨 并 且 非 常 严峻 的 现实 问 
题 。 因 而 ,研究 玻 勒 河流 域 NPP 时 空 变化 ,进而 理解 流域 的 碳 循环 过 程 和 物质 的 流动 方向 ,对 于 解决 流域 生 


http://www.ecologica.cn 


201704.00147V1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1890 生 态 学 报 37 卷 


态 环境 问题 具有 重要 的 现实 意义 。 此 外 ,研究 发 现 ,全 球 未 知 碳 汇 可 能 存在 于 陆地 生物 圈 ,特别 是 北半球 中 高 
纬度 陆地 生态 系统 ” ,而 玻 勒 河流 域 恰 好 处 于 这 一 范围 内 ,因此 对 玻 勒 河流 域 的 植被 NPP 研究 具有 重要 价 
值 。 目 前 针对 玻 勤 河流 域 长 时 间 序 列 的 NPP 估算 鲜 见 报道 , 疏 勒 河流 域 退耕 还 林 ( 草 ) .生态 综合 治理 工程 的 
实施 ,农业 灌溉 暨 移民 安置 综合 开发 项 目的 实施 ,以 及 国家 风电 基地 开发 等 引起 的 区 域 植被 NPP 空间 分 布 及 
其 年 际 变化 状况 如 何 ? NPP 对 气候 变化 可 能 的 响应 机 制 是 什么 ?这 些 问 题 或 需 回答 。 因 此 ,本 文采 用 
2001 一 2010 年 遥感 数据 .气象 站 点 观测 数据 以 及 地 面 气象 数据 ,利用 改进 的 CASA 模型 建立 玻 勒 河 流域 植被 
NPP 遥感 估算 模型 ,重建 了 流域 植被 NPP 年 际 变化 的 时 空 序 列 , 旨 在 揭示 :(1) 玻 勒 河流 域 植 被 NPP 的 空间 
分 布 格局 及 其 年 内 变化 趋势 ;(2) 气 候 因子 与 植被 NPP 间 的 关系 。 希望 本 研究 能 为 CASA 模型 在 流域 尺度 上 
的 应 用 提供 科学 依据 ,同时 为 区 域 生 态 系统 保护 措施 的 制定 和 实施 提供 决策 参考 。 


1 研究 区 与 数据 


1.1 研究 区 概况 

玻 勒 河流 域 是 甘肃 河西 地 区 三 大 内 陆 河 之 一 ,地 处 甘肃 省 河西 走廊 西 段 , 介 于 92°11’ 一 98°30'E,38? 
00' 一 42°48'N 之 间 ( 图 1) 。 行 政 区 划 上 包含 甘肃 省 的 玉 门 敦煌、 瓜 州 肃 北 县 的 大 部 分 、 阿 克 塞 县 的 一 部 分 ， 
以 及 张掖 市 肃 南 县 和 青海 省 天 峻 县 的 一 部 分 ,面积 13.36x10' km?。 干流 全 长 约 670 km, 多 年 平均 径流 量 
10.31x10* m 。 地 貌 上 由 南部 的 祁连山 地 .中 部 的 冲积 平原 和 北部 的 马 紧 山 3 大 单元 构成 , 属 典型 的 大 陆 性 若 
漠 气 候 , 年 均 气 温 6.9 一 8.8% ,大 部 分 地 区 年 均 降 雨量 在 60 mm 以 下 ,年 蒸发 量 却 高 达 3000 mm。 玻 勒 河流 域 
自 汉 代 以 来 就 一 直 是 河西 著名 的 “粮仓 ” ,灌溉 农业 发 达 。 始 于 1996 年 历经 10a 的 玻 勒 河 农业 灌溉 暨 移民 安 
置 综 合 开发 项 目 安置 了 甘肃 中 .南部 共计 7.53 万 移民 ,使 得 流域 的 人 口 和 耕地 面积 激增 ,用 水 量 加 大 ,生态 压 
力 增 大 。 由 于 玻 勒 河流 域 在 青海 省 境内 分 布 均 为 高 山 冰 川 积 雪 ,本 文 的 研究 区 不 包括 玻 勒 河上 游 在 青海 境内 
的 部 分 。 


-一 野马 > 河 
肃 北 权 大 顽 自治 县 | 
六 } 
_“ 阿 克 塞 哈萨克 族 自治 县 ~ 
Se 


92° 93° 94° 95° 96° 97° 98° 99°E 


图 1 研究 区 位 置 示意 图 


Fig.l Location of Shule river basin 


1.2 数据 来 源 
玻 勒 河流 域 面 积 较 大 , 若 采用 Landsat TM 等 中 高 分 辨 率 的 地 球 资源 卫星 数据 , 显然 不 适合 大 范围 宏观 
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研究 , 旦 数据 的 连续 性 不 强 ,时间 分 辨 率 不 高 ,难以 获得 同一 时 间 的 多 景 影像 ,影像 镶 骨 等 处 理 过 程 也 会 带 来 
一 定 的 误差 。 考 虑 到 玻 勒 河流 域 土地 覆盖 类 型 简单 ,本 文采 用 MODIS 的 归 一 化 差异 植被 指数 (NDVI) 产品 ， 
完全 能 够 满足 精度 要 求 ,而且 时 间 分 辨 率 高 ,便于 分 析 植 被 动态 变化 特征 。 数 据 获 取 自 Level 1 and 
Atmosphere Archive and Distribution System(LAADS) (http://ladsweb.nascom.nasa.gov/data/order.html) 网 站 ,为 
16d 合成 数据 ,空间 分 辩 率 为 1 km。 下 载 MOD17A3 NPP 产品 用 于 与 估算 值 进行 比较 。 气 象 数据 包括 月 平均 
气温 ,月 总 降水 量 和 月 总 太阳 辐射 , 取 自 中 国 气象 科学 数据 共享 服务 网 ( http://cdc.cma.gov.cn/home.do)。 为 
了 保证 空间 插值 精度 , 共 选 取 研 究 区 内 及 周边 的 常规 气象 站 56 个 ,太阳 辐射 监测 站 点 8 个 。 数 字 高 程 模型 
(DEM) 数 据 获 取 自 美国 地 质 勘 查 局 ( United States Geological Survey, USGS) 网 站 (http://glovis.usgs.gov/ ) ,分 
辩 率 为 1 km。 植 被 类 型 数据 来 源 于 中 国 西部 环境 与 生态 科学 数据 中 心 (http://westdc.wesitgis.ac.cn) ,为 便于 
分 析 ,将 指标 类 型 归并 为 针 叶 林 阔 叶 林 灌木 林地 .草原 .农田 .荒漠 ,水体 及 无 植被 地 段 8 类 。 


2 研究 方法 


2.1 NPP 估算 

采用 CASA 模型 计算 NPP( 表 1) ,该 模型 虽然 简化 了 生理 生态 基础 , 光 能 传递 与 转换 过 程 中 存在 一 定 的 
不 确定 性 ,但 是 该 模型 易于 实现 区 域 尺度 转换 , 且 模 型 中 的 参数 均 可 由 遥感 或 空间 分 析 方 法 获得 ,适用 于 区 域 
和 全 球 尺度 上 的 NPP 估算 。 区 域 实际 蒸 散 量 R(x , 刀 根 据 周 广 胜 5] 建立 的 区 域 实际 蒸 散 模型 求 取 ,区 域 潜在 
藻 散 量 E(x,t) 根 据 Boucher 等 ”提出 的 互补 关系 求 取 。 最 大 光 能 利用 率 s,. 的 取 值 因 不 同 的 植被 类 型 而 有 
所 不 同 ,本 文采 用 朱文 泉 等 55 和 Running 等 29 根据 生态 生理 过 程 模型 BIOME-BGG 对 10 种 植被 类 型 所 模拟 
的 值 : 针 叶 林 1.008 gCAMJ, 阔 叶 林 1.044 gCZMJ , 灌 丛 0.768 gCZMJ ,草原 0.608 gCZMJ ,农田 0.604 gCZMJ , 荒 
漠 .无 植被 地 段 和 水 体 等 取 全 球 月 平均 最 大 光 能 利用 率 0.389 gC/MJ。 


表 1 基于 改进 CASA 模型 的 NPP 计算!3"] 
Table 1 Calculation of the vegetation net primary production based on the CASA model 
编号 No。 计算 式 Calculation formula 参数 说 明 Parameter description 
APAR(x,i) 是 像 元 x 在 1 月 吸收 的 光合 有 效 辐射 (gC m 


1 NPP(x = APAR(x x 
(0) (st) tas 月 -1) ,s(x,t) 表 示 像 元 «在 1 月 的 光 能 利用 率 ( gC/MJ) 
a SOL(x,1) 表示 1 月 像 元 x 处 的 太阳 总 辐射 量 (MJ m-? 月 -1 ， 
人 FPAR 为 植被 层 对 太阳 入 射 光 合 有 效 辐射 的 吸收 比例 
3 FPAR(x,t) = aFPARAv +(1- ao)FPARw a 值 取 0.5 


NDVI4; 6) 是 像 元 * 在 1 月 的 归 一 化 差 值 植被 指数 ， 


NDVI,, , -NDVI 


i,min) 


4 FPARww = NHI TT NDVIw ww) 和 NDVIw wo 分 别 表示 第 i 种 植被 类 型 的 NDVI 
(i,max)™ (i,min) 时, 最/ FPAR FPAR | 二 0 1 和 
(FPRA,,. ~ FPRA,,) +FPRA i es .00 
| 0.95 
1 + NDVL. 
SR — SR nmin SR 为 简单 被 指数 ，SR i 
。 IPARs = 让 Ce A 为 简单 比值 植被 指数 ，SR(x,t) = 于 = NDVI ， 
ON (i,max) ON(i,min) i 别 [对 应 基 被 类 型 ,和 
(FPRAw ~ FPRAw,) + FPRAws She PU SR OA MDL HY S06 
95% 处 的 百 分 位 数 
Te1(x,t) 和 7。,(x,t) 分 别 为 低温 和 高 温 对 光 能 利用 率 的 胁迫 
6 E(x,t) = Tow,t) X Teo(x,t) X WX,t) X evax 作用 ,了 。(xz,D) 为 水 分 胁迫 影响 系数 ,sw 为 理想 条 件 下 植被 
的 最 大 光 能 利用 率 
7 T(x,t) = 0.8 +0.02 x T(x) -0.0005 x [下 Tu(*) 为 某 一 区 域 年 内 NDVI 值 达 到 最 高 时 的 月 平均 气温 
7Cx,b) 为 月 平均 温度 , 当 某 月 平均 温度 比 最 适 温度 Topt(x ) 
一 ! a 二 1 
ee 十 ee x J T(x,t) )]1} 高 10 抑或 低 13C 时 ,该 月 的 7,(x,1) 信 等 于 月 平均 浙 度 
x +e .3 x(— x)— — T(; 
1 exp ( opt(%) (x,t) ) ] | T(x,1) 为 最 适 温度 Topt(x) 时 7,,(x,t) 值 的 一 半 
9 W(x,t) = 0.5+0.5 x E(x,t)/E,(%,t) E(x,t) 是 区 域 实际 蒸 散 量 ,,(%,t) 是 区 域 潜 在 蒸 散 量 
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根据 56 个 气象 站 的 海拔 高 度 和 月 平均 气温 资料 ,计算 流域 各 月 平均 气温 的 垂直 递减 率 , 再 根据 气象 站 点 
海拔 高 度 ,将 观测 气温 换算 为 海拔 高 度 为 “0" 处 的 气温 值 。 利 用 普通 Kriging 插值 法 ,模拟 海拔 为 “0” 处 的 温 
度 7 ;根据 DEM 和 气温 垂直 递减 率 , 模 拟 起 伏地 形 下 的 真实 温度 7 , 。 降 水 数据 直接 采用 Kriging 插值 
对 气象 站 点 观测 数值 进行 插值 。 本 文 对 CASA 模型 中 计算 参数 的 主要 改进 如 下 :(1) 考 虑 到 南部 祁连山 区 和 
北部 马 聚 山地 形 复杂 ,者 仅 仅 依靠 气象 观测 站 空间 插值 的 结果 来 摘 述 流域 太阳 总 辐射 的 空间 分 布 ,误差 较 大 。 
因此 ,本 文通 过 DEM 提取 参数 ,充分 考虑 地 形 之 间 相 互 遮蔽 对 太阳 总 辐射 的 影响 ,利用 日 照 百 分 率 与 天 文 辐 
射 拟 合 太阳 总 辐射 ,获取 经 验 系数 ,从 而 较 准 确 地 计算 出 起 伏地 形 下 的 太阳 总 辐射 。(2) 本 文 在 计算 水 分 胁 
迫 因子 时 ,保持 了 原 模型 中 的 植物 生理 生态 学 因素 ,但 回避 了 其 复杂 的 土壤 参数 计算 ,只 利用 气象 数据 通过 区 
域 获 散 模 型 进行 估算。 

2.2 ”NPP 变化 分 析 方 法 
2.2.1 年 际 变化 率 计算 (趋势 分 析 ) 

一 元 线性 回归 分 析 可 以 模拟 每 个 栅 格 的 变化 趋势 ,通过 不 同时 期 单个 像 元 的 空间 变化 特征 ,综合 表征 一 
定时 间 序 列 的 区 域 格局 演变 规律 。 该 方法 的 优点 是 可 以 消除 由 于 异常 因素 对 NPP 的 趋势 分 析 所 造成 的 影 
响 , 真 实地 反映 长 时 间 序 列 NPP 的 演化 趋势 。 计 算式 为 (5 ， 


n 


n x > (i x NPP,) - Dix 和 NPP, 
Slope = i=1 i=1 i=1 (1) 


nx > 和 (> 让 
式 中 ,Slope 是 线性 拟 合 方程 的 斜率 ;NPP, 是 第 i 年 通过 最 大 值 合成 法 得 到 的 NPP 值 ,n 为 研究 时 段 的 长 度 。 
2.2.2 稳定 性 分 析 ( 变异 系数 ) 
采用 稳定 性 分 析 来 计算 2001 一 2010 年 每 个 栅 格 NPP 的 变异 系数 ,用 其 来 评估 时 间 序 列 上 的 稳定 性 , 计 


算式 为 “|. 


CG, (2) 


式 中 ,C, 为 变异 系数 ,x 为 平均 值 ,o 为 标准 差 。C, 值 越 小 ,表示 数据 分 布 越 集中 ,数据 随时 间 序 列 波动 小 , 具 
有 较 好 的 稳定 性 ;反之 则 表示 数据 分 布 较 离 散 , 数 据 随时 间 序 列 波动 较 大 。 
2.2.3 ”Hurst 指数 

采用 R/S 分 析 法 来 估算 Hurst 指数 ,从 而 分 析 NPP 的 持续 性 特征 。 篇 幅 所 限 ,计算 原理 详 见 参 考 文 
献 .。Hurst 的 值 域 为 0 一 1 之 间 , 车 0.5<H<1, 说 明 研 究 的 时 间 序 列 是 一 个 持续 性 的 序列 ,并 且 具 有 长 期 相 
关 的 特征 , 值 越 接近 于 1 ,持续 性 越 强 ; 蔡 五 =0.5, 则 表明 NPP 时 间 序 列 不 是 持续 性 序列 ,而 是 随机 性 序列 ,而 
且 它 们 之 间 并 不 存在 长 期 相关 的 特征 ; 若 0<8<0.5, 则 表现 出 反 持 续 性 ,并 且 值 越 接 近 0, 其 反 持续 性 越 强 。 
2.2.4 ”NPP 与 气候 因子 的 相关 性 分 析 

利用 基于 像 元 的 空间 分 析 法 ,计算 NPP 与 降水 或 气温 的 相关 系数 ,通过 计算 两 者 的 简单 相关 系数 ,进而 
获得 偏 相关 系数 。 计 算式 详 见 参考 文献 ”] 。 采 用 上 检验 法 对 偏 相 关系 数 进行 检验 。 


3 结果 与 分 析 


3.1 精度 验证 

在 ArcGIS 软件 支持 下 ,将 空间 化 后 的 各 参数 代入 CASA 模型 进行 栅 格 地 图 代数 运算 ,得 到 2001 一 2010 
年 逐 月 NPP 空间 分 布 。 全 球 或 区 域 尺度 上 的 NPP 佑 算 结 果 的 精度 验证 问题 ,始终 是 NPP 遥感 定量 反 演 的 难 
题 ,也 是 最 具 争 议 的 问题 。 精 度 验证 一 般 有 两 种 方法 :一 是 将 估算 值 与 实际 测量 的 数值 进行 对 比 ;二 是 将 估算 
的 NPP 空间 分 布 与 他 人 在 相同 区 域 采 用 其 他 模型 计算 的 结果 或 遥感 数据 产品 进行 对 比 。 首 先 利 用 2010 年 7 
月 和 9 月 两 次 在 玖 勒 河流 域 玉 门 . 瓜 州 . 敦 煌 和 阿 克 塞 选取 的 20 个 点 作为 野外 精度 验证 点 ,设置 样 方 为 100 
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mx100 m, 利 用 对 角 线 采样 法 :2 获得 植被 的 地 上 生产 
力 和 地 下 生物 量 ,计算 地 面 实测 NPP。 受 条 件 限制 , 样 
方 实测 地 点 多 选择 设置 在 绿洲 边缘 的 荒漠 草原 和 草 铭 
草原 内 ,农田 无 法 实测 生物 量 ,只 记录 作物 类 型 进行 转 
换 。 由 于 地 下 生物 量 实测 难度 较 大 ,部 分 样 地 的 植被 地 
下 生物 量 数据 通过 计算 根 冠 比 得 到 。 用 GPS 记录 
实测 数据 的 地 理 位 置 坐标 , 提取 相应 的 基于 改进 的 
CASA 模型 计算 出 的 NPP 数据 ,进行 佑 算 值 的 精度 验 
证 。 相 关 性 分 析 ( 图 2) 表 明 ,估算 值 与 地 面 实测 值 之 间 
的 相关 性 非常 显著 ,决定 系数 R* =0.8811,P<0.01, 说 明 
改进 的 CASA 模型 适用 于 踊 勒 河流 域 NPP 的 估算 。 误 
差分 析 显 示 ,CASA 模型 估算 的 NPP 与 地 面 实测 NPP 
的 平均 绝对 误差 为 37.13 gC m” a ' ,平均 相对 误差 为 
0.208 , 均 较 小 。 因 实际 样 方 大 小 (100 m) 与 MODIS 像 
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图 2 NPP 实测 值 与 模型 估算 结果 的 比较 


Fig.2 Comparison between measured and estimated NPP 


元 分 辨 率 (1 km) 差距 很 大 ,误差 在 所 难免 。 总 体 来 看 , 佑 算 值 在 荒漠 草原 偏 低 而 在 草 多 草原 偏 高 。 


虽然 野外 调查 所 获取 的 生物 量 数 据 比 较 可 靠 , 但 很 难 在 整个 流域 区 内 进行 大 范围 相对 均匀 的 实地 调查 取 
样 。 加 之 研究 时 段 跨 度 为 10a, 生 物 量 分 布 的 空间 和 时 间 异 质 性 较 大 ,因此 如 果 简 单 地 利用 有 限 的 实地 调查 
所 获得 的 平均 生物 量 来 推算 整个 区 域 的 生物 量 , 则 可 能 产生 较 大 误差 。 为 此 ,本 文 将 估算 结果 与 2001 一 2010 
年 MODIS 的 1 km 年 分 辨 率 NPP 产品 (MOD17A3 ) 进行 了 对 比 ,发 现 估算 的 NPP 与 MOD17A3 产品 在 空间 分 
布 上 较为 吻合 ,除去 2001.2003 .2007.2010 年 这 4a 中 MOD17A3 值 与 估算 值 相差 稍 大 外 ,其余 年份 相差 很 小 ， 
差 值 介 于 2.18 一 21.53 gC m ”a 之 间 。MOD17A3 年 总 量 与 本 文 估算 结果 总 体 差距 亦 不 大 , 差距 较 大 的 年 份 
是 2001 年 和 2003 年 , 且 除 2001 年 外 ,各 年 份 估算 的 NPP 总 量 均 略 高 于 MOD17A3 值 。 总 体 而 言 ,MOD17A3 


的 NPP 年 际 波动 明显 大 于 本 文 估算 值 。 
3.2 年 平均 NPP 的 空间 分 布 

2001 一 2010 年 , 玻 勒 河流 域 年 均 单 位 面积 上 的 
NPP 空间 分 布 状况 见 图 3。 玻 勒 河流 域 NPP 值 小 于 
100 gC ma 的 区 域 占 流域 总 面积 的 比重 高 达 62.9% ， 
反映 出 流域 整体 生产 力 水 平 不 高 。 从 空间 分 布 来 看 , 流 
域 NPP 值 呈 现 南 高 北 低 、 东 高 西 低 的 态势 ,表现 出 较为 
明显 的 经 向 地 带 性 “条 带 ” 分 布 特征 , 即 上 游 的 祁连山 
区 和 中 游 的 绿洲 区 NPP 较 高 ,其 两 侧 的 戈壁 荒漠 和 洪 
冲积 扇 NPP 较 低 。 这 也 是 干旱 内 陆 河流 域 普 遍 具 有 的 
特点 ,体现 了 绿洲 一 荒漠 一 山地 耦合 生态 系统 的 典型 特 
征 。 流 域 NPP 最 高 的 区 域 是 上 游 地 带 , 西 祁连山 高 寒 
荒漠 草原 区 . 柴 达 木 盆地 东北 部 山地 高 寒 荒漠 草原 区 以 
及 阿尔 金山 高 寒 荒漠 草原 区 的 NPP 值 一 般 均 在 150 gC 
ma 以 上 ,自然 植被 以 高 山 草 甸 为 主 。 流 域 中 下 游 地 
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3 ”2001 一 2010 年 琉 勒 河流 域 NPP 平均 值 的 空间 分 布 
Fig.3 Spatial distribution of 10-year mean annual NPP of Shule 
River basin in 2001 一 2010 


势 低 平 ,NPP 高 值 区 均 分 布 在 以 灌溉 农业 为 主 的 绿洲 区 。 


北部 的 马 甘 山 风 蚀 荒 漠 区 和 花 牛 山 - 柳 园 淮 木 - 半 灌 


木 蕊 漠 区 植被 稀 疏 ,大 部 分 区 域 NPP 在 50 一 100 gC m” a 之 间 。NPP 最 低 的 区 域 分 布 在 敦煌 西部 的 库 木 塔 
格 沙漠 . 阿 克 塞 北部 的 安南 坝 以 及 瓜 州 的 安西 极 旱 荒漠 保护 区 ,NPP 值 大 都 在 50 gC ma 以下。 
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3.3 NPP 年 际 变化 特征 

玻 勤 河流 域 年 均 NPP 值 2001 年 最 低 ,2010 年 最 NPP 什 /gC m2 a 
高 ,分 别 为 99.97 gC m- ar 和 127.2 gC ma ',2001 一 人 
2010 年 平均 NPP 值 为 125.08 gC m ”a ”"。 从 NPP 值 的 100 1150 ~ 
变化 来 看 (图 4) ,大 致 可 分 为 两 个 阶段 : 2001 一 2008 。 、 go 0 
年 ,NPP 由 急剧 上 升 变 为 2008 年 之 后 的 小 幅 波 动 ， 妾 。 | 二 | Ee 
2008 一 2010 年 ,NPP 由 急剧 下 降 转 为 急剧 上 升 。 整 体 RE | | | | | / ie 
来 看 ,10 间 NPP 值 整体 呈现 出 波动 中 上 升 的 趋势 。 从 于 8 lmse 
2001 一 2010 年 ,NPP<50 gC ma 的 区 域 面积 整体 上 20 | | | | | | | | | 100 
呈 减 少 的 趋势 ,而 NPP>350 gC m? a 的 区 域 面积 整体 os 
则 呈现 增加 的 趋势 。 8 

年 份 Year 


3.4 “NPP 空间 变化 特征 
3.4.1 总 体 趋势 

图 5 为 10 年 来 琉 款 河 流域 年 最 大 NPP 的 变化 趋 
势 ,流域 年 最 大 NPP 呈 缓 慢 增加 趋势 ,平均 趋势 系数 为 
1.79。NPP 增加 的 区 域 占 总 面积 25.15% ,其 中 显著 增 


图 4 2001 一 2010 年 疏 勒 河流 域 年 平均 NPP 变化 趋势 
Fig. 4 
Northwest Arid Region from 2001 to 2012 


The change trend of monthly average NPP in the 


加 区 域 为 6.59% , 轻 度 增加 区 域 为 18.56% , 主要 分 布 于 南部 的 肃 北 、 玉 门 的 南部 及 阿 克 塞 的 东部 ;基本 不 变 的 
区 域 约 占 62.92% ;减少 的 区 域 较 少 , 约 占 总 面积 的 11.93% ,其 中 轻 度 减少 区 域 为 10.11% ,显著 减少 的 区 域 仅 


为 1.81%, 主 要 分 布 于 北部 的 肃 北 、 阿 克 塞 中 部 及 玉 门 、 瓜 州 .敦煌 的 绿洲 周围 。 
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5 2001 一 2010 年 疏 勒 河流 域 NPP 线性 变化 趋势 图 
Fig.5 ” Distribution of linear trend of annual NPP changes in the Shule River basin from 2001 to 2010 


3.4.2 NPP 变化 的 稳定 性 
由 10 年 间 下 b 勒 河流 域 NPP 变异 系数 分 布 图 (图 6) 


可 知 , 玖 勒 河流 域 NPP 整体 呈现 


头 大 一 头 小 ,高 低 


波动 并 存 , 低 度 波动 居多 ”的 分 布 格局 。 整 个 流域 不 同 地 区 的 NPP 稳定 性 存在 明显 的 空间 差异 ,高 波动 区 多 


分 布 在 白杨 河 石油 河 . 玻 勒 河和 干流、 榆林 河 .党 河 .安南 坝 河 等 河流 两 侧 , 说 明 这 些 区 域 易 受 到 自然 或 人 为 扰 
动 , 植 被 的 稳定 程度 低 。 中 等 波动 区 域 大 都 分 布 在 上 游 和 中 下 游 之 间 的 冲 、 洪 积 扇 。 低 波动 区 域 则 连 片 分 布 
于 流域 北部 和 西部 。 低 波动 区 域 占 流域 总 面积 的 56.74% , 较 低 波动 区 占 38.12% , 而 高 波动 区 仅 占 0.3% ,这 


说 明 下 勤 河流 域 NPP 整体 相对 稳定 。 


http://www.ecologica.cn 


201704.00147V1 


chinaXiv 


ChinaXiv 合 作 期 刊 


6 期 潘 竟 虎 ”等 .2001 一 2010 年 玖 勒 河 流域 植被 净 初 级 生产 力 时 空 变化 及 其 与 气候 因子 的 关系 1895 


3.4.3 ”NPP 变化 的 未 来 趋势 特征 

前 述 分 析 主 要 是 对 10 年 来 玻 勒 河流 域 NPP“ 格 局 - 
过 程 "的 解释 ,而 对 于 其 未 来 的 变化 趋势 尚 不 明确 ,为 
此 ,借助 于 Matlab 和 ArcGIS 软件 ,计算 研究 区 NPP 的 
Hurst 指数 。Hurst 的 值 域 范 围 为 0.1 一 1, 均 值 为 0.58。 
反 持 续 性 占 流 域 总 面积 的 21.70% ,持续 性 占 流 域 面 积 
的 78.30% ,说 明 玻 勒 河流 域 NPP 变化 的 同 向 特征 要 强 
于 反 向 特征 。 将 Hurst 指数 划分 为 弱 、 中 、 强 3 个 持续 
性 类 型 ,其 阔 值 分 别 为 :<0.5,0.5 一 0.65 和 >0.65。 从 
弱 .中 、 强 的 程度 来 看 ,中 强 程度 的 反 持 续 性 仅 占 流域 面 91° 92°% 9” 中 9 96°% 97 98 99°E 
积 的 3%, 弱 反 持 续 性 占 18.7%, 弱 持续 性 占 49.93% ,说 图 6 ”2001 一 2012 年 疏 勒 河流 域 NPP 变化 稳定 程度 
明 朴 勤 河 流域 有 一 半 的 区 域 未 来 NPP 的 变化 与 过 去 Fig.6 Standard deviation of NPP changes in the Northwest Arid 
NPP 的 变化 相 一 致 ,但 这 种 持续 性 表现 得 并 不 够 明显 。 Region from 2001 to 2012 
中 、 强 持续 性 所 占 比重 分 别 为 19.83% 和 8.549% , 这 些 区 
域 未 来 NPP 的 变化 与 过 去 NPP 变化 有 较 强 的 相关 性 。 
3.4.4 ”NPP 变化 的 趋势 性 和 持续 性 综合 分 析 

将 一 元 线性 拟 合 的 Slope 值 与 Hurst 指数 计算 结 
进行 琶 加 分 析 , 得 到 2001 一 2010 年 NPP 变化 趋势 的 持 
续 性 分 布 (图 7) ,为 使 结果 更 加 清晰 明了 , 共 分 为 12 种 
情形 。 玻 勒 河流 域 不 变 与 弱 持 续 性 .不 变 与 强 持续 性 、 
增加 与 弱 持续 性 ,不 变 与 反 弱 持续 性 占 研 究 区 面积 的 百 
分 比分 别 为 30.50% .18.02% .12.54% 和 12.15% ,可 见 整 
个 流域 在 10 年 间 NPP 分 布 并 未 大 幅度 地 发 生变 化 , 且 
可 持续 性 较 差 ,这 与 前 文 分 析 的 下 勒 河流 域 NPP 变化 
趋势 一 致 。 综 合 分 析 ,流域 减少 且 未 来 会 持续 性 减少 的 
面积 为 9.15% , 主要 集中 在 肃 北 . 阿 克 塞 . 教 煌 地 区 ,应 9 
加 强 这 些 地 区 的 生态 环境 治理 ;增加 旦 未 来 会 继续 增加 号 11 城 少 与 强 反 持续 性 。” 国 23 不 变 与 弱 持续 性 

国 12 减 少 与 轮 反 持续 性 。 国 24 不 变 与 强 持续 性 


的 面积 为 20.63% ,主要 分 布 于 南部 肃 北 、 玉 门 瓜 州 的 南 国 13 减 少 与 弱 持 续 性 国 31 增 加 与 强 反 持续 性 
国 14 减少 与 强 持续 性 0 32 增 加 与 弱 反 持续 性 


93° 94?” 95° 96° 97? 98"” 99° 100°E 


部 和 阿 克 塞 的 东部 ,这 些 地 区 基本 位 于 玻 勤 河 流域 的 上 国 | 21 不 变 与 强 反 持 续 性 ”上 轩 33 增 加 与 弱 持续 性 
游 ;特别 应 注意 的 是 不 变 且 为 弱 持 续 性 的 地 区 ,其 比例 0 


高 达 30.50%, 由 于 可 持续 性 较 差 的 原因 ,应 防止 其 向 退 图 7 ”基于 趋势 性 和 Hurst 指数 的 NPP 空间 变化 特征 图 
化 的 趋势 发 展 ,这 些 地 区 主要 分 布 于 流域 的 西北 部 。 可 Fig.7 The spatial change characteristics of NPP based on the 
以 预见 玻 勒 河流 域 NPP 值 在 未 来 的 几 年 还 是 会 以 不 变 Slope and Hurst index 
为 主 ,增加 为 辅 ,但 不 变 的 可 持续 性 较 差 , 应 继续 加 强生 
态 环 境 的 治理 。 
3.5 ”NPP 与 气候 因素 的 相关 关系 
3.5.1 年 NPP 与 年 气候 因子 的 相关 关系 

植被 的 生长 与 气候 因素 密切 相关 ,由 图 8 可 知 , 琉 勒 河流 域 植被 NPP 与 年 降水 和 年 均 温 的 平均 相关 系数 
分 别 为 0.37 和 0.02 ,平均 偏 相 关系 数 分 别 为 0.41 和 -0.05 ,说 明 踢 勤 河流 域 植 被 年 NPP 总 体 上 与 年 降水 的 相 
关 性 更 高 ,呈正 相关 关系 。 比 较 而 言 ,流域 西北 部 受气 温 的 影响 较 大 ,而 东南 部 受降 水 的 影响 大 。 图 8 是 NPP 
与 降水 的 相关 系数 分 布 图 ,消除 气温 影响 后 (图 8) ,研究 区 内 NPP 与 降水 呈现 正 相关 性 的 区 域 面积 增 大 , 且 
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相关 性 增强 ; 而 呈现 负 相关 性 的 区 域 面 积 减 小 ,相关 性 减弱 。NPP 与 气温 的 相关 系数 空间 分 布 如 图 8 所 示 ， 
消除 降水 影响 后 (图 8) ,研究 区 内 NPP 与 气温 呈现 正 相关 性 的 区 域 面积 减 小 ,而 呈现 负 相关 性 的 区 域 面 积 增 
大 ,整体 相关 性 增强 。 将 相关 系数 与 偏 相 关系 数 空间 分 布 图 对 比 ,发 现 大 部 分 地 区 的 偏 相 关系 数 较 之 相关 系 
数 有 所 下 降 ,说 明 这 些 地 区 受降 水 和 气温 的 共同 影响 。 


8 年 NPP 与 年 降水 和 年 均 温 的 相关 系数 及 偏 相 关系 数 


Fig.8 Correlation relationship between annual NPP and annual precipitation and annual mean temperature 


3.5.2 月 NPP 与 当月 气候 因子 的 相关 关系 

玻 勤 河流 域 月 均 NPP 与 当月 降水 及 当月 均 温 的 平均 相关 系数 分 别 为 0.81 和 0.93 ,平均 偏 相 关系 数 分 别 
为 -0.13 和 0.81, 可 见 月 NPP 与 当月 气候 因子 的 显著 性 都 较 年 NPP 与 气候 因子 的 相关 性 要 高 。 肃 北 县 的 南 
部 受降 水 的 影响 较 大 ,其 他 县 区 主要 受气 温 的 影响 较 大 。 月 NPP 与 当月 降水 量 成 正 相 关 的 面积 占 研 究 区 面 
积 的 34.98%( 图 9) ,月 NPP 与 当月 气温 成 正 相 关 的 面积 占 人 研究 区 面积 的 99.68% ,说 明 植 被 NPP 在 月 尺度 变 
化 上 ,温度 是 其 主要 影响 因子 。 
3.5.3 月 NPP 与 前 一 月 气候 因子 的 相关 关系 

气候 因子 对 植被 的 生长 具有 一 定 的 时 滞 效 应 ,对 于 不 同 地 区 .不同 的 气候 因子 时 清 效 应 程度 也 不 同 。 由 
图 10 可 知 ,流域 月 NPP 与 前 一 月 降水 量 、 前 一 月 的 月 均 温 的 平均 相关 系数 分 别 为 0.67 和 0.69 ,平均 偏 相 关系 
数 分 别 为 0.20 和 0.27, 可 见 月 NPP 与 前 一 月 的 降水 和 当月 降水 的 相关 系数 的 大 小 接近 ,而 月 NPP 与 前 一 月 
的 气温 的 相关 系数 较 当 月 的 气温 相关 系数 有 明显 的 下 降 。 对 整个 距 勒 河流 域 而 言 , 月 NPP 对 气候 因子 变化 
的 时 滞 效 应 并 不 明显 ,但 从 空间 分 布 来 看 , 肃 北 南部 敦煌 及 阿 克 塞 的 部 分 地 区 的 NPP 对 降水 的 时 滞 效 应 较 
为 明显 ,而 在 肃 北 北部 和 玉 门 ,NPP 对 气温 的 时 滞 效 应 更 为 明显 。 
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9 月 NPP 与 当月 降水 和 均 温 的 相关 系数 及 偏 相 关系 数 


Fig.9 Correlation relationship between monthly NPP and monthly precipitation and monthly mean temperature 


4 讨论 
目前 ,前 人 关于 琉 勒 河流 域 NPP 的 模拟 研究 极 少 , 而 且 模 拟 的 疏 惑 河流 域 不 同 生态 系统 NPP 值 存在 很 


大 差异 。 本 文 模拟 的 2001 一 2010 年 疏 款 河 流域 NPP 平均 值 为 102.26 gC ma ,与 张 娟 等 5 的 结果 (95.2 
gC ma !) 接 近 , 高 于 刘 春 雨 等 .的 模拟 结果 (55.89 gC ma:) 。 就 植被 类 型 而 言 ,本 文 估算 的 安西 荒漠 草 
地 NPP 平均 值 为 266.47 gC m? a ,也 与 张 娟 等 "的 和 张杰 等 :55 的 研究 结果 接近 。 不 同学 者 模拟 结果 不 一 臻 
的 主要 原因 :1) 研究 的 时 间 序 列 不 一 致 ;2) 数 据 源 的 时 空 分 辩 率 不 同 ;3) 所 用 的 模型 不 同 ;4) 使 用 CASA 模型 
时 ,气象 数据 的 插值 大 都 未 考虑 地 形 对 气温 .降水 和 太阳 辐射 空间 分 布 的 影响 。 

2001 一 2010 年 玻 勒 河流 域 交 初级 生产 力 的 空间 变化 以 波动 式 增 加 为 特征 ,NPP 增加 显著 且 持 续 性 较 强 
的 区 域 主要 分 布 在 肃 北 县 南片 区 阿 克 塞 县 的 玖 勒 河 上 游 地 区 ,以 及 中 游 的 玉 门 和 瓜 州 绿洲 农田 区 ;而 NPP 
显著 减少 上 且 持续 性 较 强 的 区 域 主要 分 布 于 肃 北 北 片区 的 马 聚 山地 区 . 玻 勒 河 下 游 干 流 以 及 各 支流 的 绿洲 外 于 
地 区 。 在 干旱 内 陆 河流 域 ,水 资源 是 维系 绿洲 一 荒漠 生态 系统 发 育 与 演化 的 首要 因素 。 玻 勒 河中 下 游 绿洲 地 
区 农田 为 主 的 NPP 持续 增加 与 绿洲 外 围 荡 漠 草 地 植被 NPP 的 持续 减少 ,在 一 定 程 度 上 反映 出 人 类 活动 对 水 
资源 的 调控 在 农业 生产 上 的 积极 效应 ,但 另 一 方面 ,农业 用 水 .生产 用 水 .生活 用 水 的 增加 使 得 生态 用 水 被 严 
重 挤占 ,引起 自然 植被 退化 和 生产 力 下 降 。 仅 玻 勒 河 下 游 的 敦煌 市 和 瓜 州 县 ,人 口 就 由 1989 年 的 18.26 万 增 
加 到 2010 年 的 33.48 万 。 尽 管 气候 因素 是 疏 勒 河流 域 NPP 时 空 变化 的 重要 影响 因素 ,但 研究 时 段 内 人 类 活 
动 的 不 断 增强 对 植被 NPP 的 影响 也 不 断 加 大 ,长远 来 看 , 必 将 反作用 于 人 类 生活 ,产生 不 利 影响 。 按 照 “ 反 规 
划 ” 的 思想 ,从 流域 NPP 变化 及 其 影响 因子 的 空间 分 布 情况 ,可 以 判断 出 某 区 域 影响 植被 NPP 演变 的 限制 因 
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图 10 月 NPP 与 前 一 月 降水 和 均 温 的 相关 系数 及 偏 相关 系数 


Fig.10 ”Correlation relationship between monthly NPP and previous monthly precipitation and monthly mean temperature 


素 ,从 而 对 该 区 采取 有 目的 的 人 为 干预 ,确保 植被 恢复 与 生态 重建 ,维护 区 域 生 态 安 全 。 

NPP 的 时 空 变化 受 多 种 因素 的 综合 影响 ,本文 所 采用 的 改进 CASA 模型 的 模拟 结果 取决 于 植被 指数 和 气 
象 数据 ,未 考虑 土壤 状况 ,大气 二 氧化 碳 浓度 以 及 人 类 活动 等 对 NPP 的 影响 ,计算 结果 与 实际 状况 必然 存在 
偏差 ,但 遥感 反 演 的 NPP 是 一 个 综合 性 指数 ,可 以 反映 植被 长 势 与 生境 条 件 , 故 本 文 的 研究 结果 仍 具有 科学 
参考 价值 。 目 前 已 经 提出 的 NPP 估算 模型 和 方法 ,各 有 其 优势 和 不 足 , 尚 无 统一 的 普 适 性 方法 。 本 文 虽 对 
CASA 模型 进行 了 改进 ,但 所 采用 模型 中 有 些 静 态 参 数 是 直接 利用 已 有 的 研究 成 果 , 这 些 参数 多 是 依赖 于 较 
小 尺度 的 实验 区 或 少数 试验 站 的 监测 数据 测算 得 到 的 ,直接 将 其 利用 到 大 区 域 上 ,必然 会 对 佑 算 结果 产生 影 
响 。 对 于 各 种 类 型 的 区 域 . 不 同 的 尺度 采用 适合 建 模 参数 将 是 下 一 步 研 究 的 重点 。 此 外 ,研究 区 气象 站 点 
(尤其 是 太阳 辐射 站 点 ) 较 少 且 分 布 不 均 ,虽然 对 气象 数据 插值 方法 做 了 改进 ,但 是 还 需 发 展 更 加 有 效 和 精确 
的 插值 方法 以 提高 模型 估算 精度 。 
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